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We present evidence for unquenching effects in Nf = 2, 16'^32 ensembles by comparing with 'equivalent' quenched 
data at fo — 5.0. A (small) VEV for torelons signals (weak) string breaking. A 15 — 20% reduction in the scalar 
glueball mass relative to quenched is argued to be (in part at least) a discretisation effect. We find a chiral 
suppression of the topological susceptibility consistent with expectations, and agreement between fermionic and 
gluonic methods for measuring the topological charge. 



To study the effects of vacuum polarisation in 
QCD, we should vary the sea quark mass, m^, 
having first removed discretisation effects. With 
confident continuum extrapolations of many lat- 
tice quantities unavailable, hovi^ever, we must in- 
stead try to 'freeze' lattice artefacts. UKQCD 
does this by attempting to fix the lattice Som- 
mer scale, fo ~ 5, as to, is changed. Variations 
due to residual changes in f o are further reduced 
by the the leading order discretisation errors for 
the C'(a)-improved action being quadratic, rather 
than linear, in the lattice spacing, a. An 'equiva- 
lent' value of fo is found in the quenched, 'gluody- 
namical' theory at coupling /3wii = 5.93. Details 
of scales, parameters and matching can be found 
in |l]. We use as a measure of rriq the square of 
the pseudoscalar 'pion' mass. Circumflexes de- 
note dimensionless lattice quantities. 

The failure of experiments to detect unam- 
biguous glueball states is believed to be a conse- 
quence of mixing between the light glueballs and 
qq states. A major goal of lattice QCD is to pre- 
dict the mixing matrix elements, but the uncer- 
tainties in attempts so far |^,|| mean, however, 
that phenomenological attempts to describe the 
content (gluonic or qq) of the scalar sector glue- 
ball candidates lead to widely differing results [ 4 1 . 

As a preliminary to a full mixing analysis [5|, 
we study the light glueballs and winding flux 
tubes on UKQCD dynamical configurations and 
find, in common with , that accurate mass esti- 
mates are possible on moderately sized ensembles 
(4000 - 8000 HMC trajectories). Full results can 



be found in @] , along with a discussion of the ap- 
plicability of the term 'glueball' when mixing is 
induced by vacuum polarisation. 

On a spatial torus, Polyakov loops (PL) couple 
to colour flux tubes that close upon themselves 
through the periodic boundary with unit winding 
number. In the conflning phase of gluodynamics 
such 'torelons' are stable, with the PL vacuum 
expectation value (VEV) zero even after operator 
improvement. To a first approximation the flux 
tube mass is the product of the string tension, K, 
and the spatial extent of the lattice, L. 

Dynamical quarks break the centre symme- 
try that ensures the orthogonality of contractible 
and non-contractible operators. The PL can 
now have a non-zero VEV. This VEV should in- 
crease as rriq — > (and centre symmetry break- 
ing gets worse) and, due to the spatially periodic 
fermionic boundary conditions, be negative. In 
a simple model it goes as — (m,)^-^ for (moder- 
ately) heavy sea quarks This is 'string break- 
ing' of the confining flux tube, and Fig. |^ confirms 
this picture for improved PL VEVs. With the 
breaking of centre symmetry, mixings between 
torelons winding in different directions becomes 
possible. Also allowed is mixing between single 
(non-contractible) PL and (contractible) glueball 
operators. This should be the leading source of 
(mass reducing) flnite volume effects in the scalar 
glueball: the double torelon that plays this role in 
the quenched theory is too heavy to figure here. 
Both mixings are very weak, however, as expected 
for a numerically small VEV, and we can be con- 
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Figure 1. The VEV of the spatial torelon. 

fident that finite volume effects are small for the 
scalar glueball. We can also thus use the torelon 
mass to obtain K, which agrees well with that 
from static interquark potential fits ^ . 

The masses of the lightest scalar and tensor 
states are shown in Fig. ^, with the 'equiva- 
lent' quenched results at /3wii = 5.93, and the 
threshold for the scalar to decay to pseudoscalars. 
Whilst the (admittedly noisy) tensor mass shows 
no significant quenching effect, vacuum polarisa- 
tion reduces the scalar mass by 15 — 20%. The 
at most weak dependence of this figure on rrig 
and the nn threshold suggest it is an increased 
discretisation effect relative to the quenched the- 
ory, and that it may not persist in the continuum 
limit. In gluodynamics, this 'scalar dip' is (prob- 
ably) due to the scaling trajectory, /3wii — * oOi 
passing close to a critical point in the extended 
fundamental-adjoint gauge coupling plane. The 
further reduction in the presence of sea quarks 
for the C'(a)-improved action is then not unex- 
pected: hopping parameter expansions suggest 
the scaling trajectory now passes even closer to 
this critical point. Measurements of the couplings 
in such a mixed effective action appear to support 
this. The lessened scalar mass at fixed, finite a 
is thus not necessarily indicative of a continuum 
effect 01. 
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Figure 2. The lightest glueball masses. 

The ability to vary parameters fixed in Nature 
has made lattice Monte Carlo simulation a valu- 
able tool for investigating the role of topological 
excitations in QCD and related theories. 

The topological susceptibility, x = {Q^)/V: is 
the squared expectation value of the topological 
charge. 
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normalised by the volume. Sea quarks induce 
an instanton-anti-instanton attraction which be- 
comes stronger in the continuum limit, suppress- 
ing X 0]- Given various assumptions [p^, in 
the chiral limit of Nf fiavours we have x = 
f^m1/{2Nf) + 0{m^). This is in a convention 
where the experimental value of the pion decay 
constant ~ 93 MeV, and requires f^m'^V » 
1, which holds for the lattices used here. The 
topological susceptibility is measured on UKQCD 
gauge fields using a symmetrised 'Wilsonian' op- 
erator, Q{n), having first performed 10 Wilson ac- 
tion cools to regulate the ultraviolet dislocations. 
Fig. P shows the chiral suppression in x relative 
to the quenched limit with an interpolating fit ||^ . 
Including the term, we can fit the four most 
chiral data, obtaining — 105 ± 6 t-io MeV at 
a ~ 0.1 fm. A fuller discussion of these results is 
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Figure 3. The topological susceptibility. 

available in ||9|,|l0|. 

The topology of the vacuum is clearly influ- 
enced by the sea quarks. At finite a, however, 
Q is method dependent. To understand the sys- 
tematic biases this introduces, we compare sev- 
eral contrasting methods with differing system- 
atic errors for measuring Q on a QCDSF ensem- 
ble for a ~ 0.1 fm. The first is as described above. 
We also use the '5Li' method ||l3|, which is ex- 
act up to 0{a^) and approximately restores scale 
invariance for large lattice instantons. We use 
the Boulder smearing technique, which mimics a 
renormalisation group inspired smoothing of the 
gauge fields | p^ . Finally, we employ a fermionic 
charge operator on the 'hot' gauge fields p5[ |. 

Fig. ^ shows Monte Carlo time series of the var- 
ious Q, which have not been rounded to the near- 
est integer. We can quantify the visually striking 
agreement through normalised correlation func- 
tions using pairs of measurements. Preliminary 
results give 85 — 95% agreement between any pair 
of charge definitions; subsequent work will com- 
pare systematic effects and the statistical accu- 
racy of topological susceptibility measurements. 
Other work in this area can be found in [p^ . 

We thank Ph. de Forcrand for access to his 
code, and T. Kovacs and D. Pleiter for results 
used in Fig. H. 
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Figure 4. Comparing topological charges. 
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